On the applicability of the two-band model to describe transport across n-p 

junctions in bilayer graphene 

C. J. Poole 

Department of Physics, Lancaster University, Lancaster, LAI 4YB, UK 



o 

C 



Abstract 

We extend the low-energy effective two-band Hamiltonian for electrons in bilayer graphene (Ref . [T] ) to include a spatially 
dependent electrostatic potential. We find that this Hamiltonian contains additional terms, as compared to the one used 
earlier in the analysis of electronic transport in n-p junctions in bilayers (Ref. [3]). However, for potential steps |u| < 71 
(where 71 is the interlayer coupling), corrections to the transmission probability due to such terms are small. For the 



angle-dependent transmission T(8) we find T{6) 
factor: F = 0.241 for u = 40meV. 



sin (20) — (2u/37i) sin(4f?) sin(0) which slightly increases the Fano 
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Graphene, a crystal of carbon atoms in a two-dimensional 
(2D) honeycomb lattice, is a gapless semiconductor [HE]. 
Gating of graphene enables one to vary the carrier density 
and therefore move the Fermi level from the conductance 
band to the valence band. Gating graphene flakes with 
multiple gates enables one to generate electrostatically de- 
fined n-p junctions [51415]. Bilayer graphene in particular is 
often described by a four-band Hamiltonian from a tight- 
binding calculation (given that there are four atoms in the 
unit cell; see Fig. [T]). For low energies near the Fermi sur- 
face, one can describe the transport of electrons with a 
two-band Hamiltonian pQ. Transport across an n-p junc- 
tion in bilayer graphene in the low-energy, ballistic regime 
has been previously studied in Ref. [5J, but without con- 
sidering the possibility of a correction due to the spatial 
dependence of the electrostatic potential. 

In this paper, we extend the derivation of an effec- 
tive two-band Hamiltonian for bilayer graphene (in the 
low-energy regime) to include the effects of a spatially de- 
pendent electrostatic potential u, and a gap in the energy 
spectrum A. The re-derived two-band model Hamiltonian 
contains several additional terms which originate from the 
spatial derivatives of u{x). We use this in the analysis of 
the problem of an n-p junction, where we find a change 
in transmission probability, as compared to the analysis 
in Ref. [5] , which showed perfect transmission through the 
n-p junction at an angle of 45° (see Fig. [3|. This analysis 
shows that the additional terms in the effective two-band 
Hamiltonian induced by the gradient expansion involving 
the lateral potential are small, and thus the correctional 
term to the angular transmission probability increases the 
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Figure 1: Schematic of AB (Bernal) stacked bilayer graphene show- 
ing intralayer and interlayer couplings, as well as a unit cell compris- 
ing of four carbon atoms: A,B,A2,B2. Inset: energy bands in bilayer 
graphene near a Kpoint. The energy of the quasiparticles is near 
e = 0, qualifying the assumption that 71 is large compared to other 
energies in the system. The transformation reduces the band struc- 
ture to blue (solid) bands only. 



angle at which perfect transmission occurs by a few de- 
grees. This also results in a small correction to the Fano 
factor. 

Using the nearest-neighbour tight-binding approxima- 
tion in the Slonczewski-Weiss-McClure parameterisation 
[TS] , one can write the Hamiltonian at a K point (for basis 
((I>a,<I>B2,4>A2,4>b)) as 
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where cr = (a x ,a y ), p = (p x ,p y ) and £ is the Dirac point 
index (£ = +1 for the valley around the K point, —1 for 
the valley around the K' point, and throughout this paper 
we set H = 1). v = ^ajn/h and cr.; are the Pauli spin ma- 
trices. Furthermore, A = £2 — £1 is the difference between 
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Figure 2: Low-energy band structure of a single valley on either side 
of the potential step. The Fermi energy is the same on both sides, 
causing an electron from the left side to tunnel through the barrier 
from the conductance band to the valence band on the right side. 



the on-site energies in the two layers, £2 = |A, E\ — — |A, 
which produces a gap in the energy spectrum [17] . A po- 
tential term u — \u is added along the diagonal to rep- 
resent the electrostatic potential (we neglect inter-valley 
scattering between K and K'; I is the unit matrix). We 
assume that the interlayer coupling ji is large compared 
to other energies in the system (which is reasonable for 
the low-energy regime near the Dirac points). Given that 
s <C 71 (where e is the energy of charge carriers) , and with 
e = p 2 /2m, where m = -fi/2v 2 , we see that (fw/71) 2 <C 1. 
From this justification, we drop terms beyond quadratic 
in momentum in the following calculations. We assume a 
non-adiabatic system, with 



a< l ± < (Z,A F ), 



(2) 



where a is the lattice constant, I the width of the step (see 
Fig. [2]), l± — u/71, and Xp is the Fermi wavelength. 

We use a Schrieffer- Wolff transformation [T5] to map 
Eq. ([lj in a 4D Hilbcrt space into a 2D subspace, creating 
an effective Hamiltonian. If we let %4 X 4 = H° + 6H, with 
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we can then write the associated Green's function as (? 4x4 = 
(e — H4X4) = (e — H° — SH)^ 1 and expand: 



g 4x4 = (e - n )- 1 + (e - H°)~ 1 8'H(e - H^ 1 - 



(4) 



Given the basis that %4 X 4 is constructed in, and that 
the low-energy quasiparticle transport is directly from atom 
A to B2 in the bilayer unit cell [T] (see Fig. [T]), we wish to 
map "^4x4 onto the Hu block matrix, using a Schrieffer- 
Wolff transformation. This has the effect of only keeping 
terms with an even number of SH components. The end 
result is that G^ 1 = s — Hu — Hi 2 (e — H 2 2)~ 1 H2i- 

During this projection, the orthonormality of the wavevec- 
tors has to be preserved. To do this, we notice that G^ 1 
is an inverse Green's function of the form 



G u ' = e - Hu + n + e/3, 

P = —0H12H21, 
7J 
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We write an effective Schrodinger equation as e(l + 
(3)\ip) = {Hu - Cl)\ip). We wish to enforce (ip\l + /3\xJj) = 
\ip\ 2 and ((p\ip) — 0. Writing the wavefunction in terms of a 
new wavefunction = (I + /3) _1 / 2 |-0). Inserting this 

result back into the effective Schrodinger equation gives us 
H cff = (l + /3)- 1 / 2 (iJii-fi)(l + ^)- 1 / 2 , which after Taylor 
expanding around /3 up to 0{(3 2 ) produces 



H eS = Hu-tl-{(Hu-n),^H 12 H 21 }, (6) 
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where the curly braces denote the anticommutator. The 
effective Hamiltonian can thus be calculated as 
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The first two lines form the Hamiltonian found in Rcf. 
P] (neglecting trigonal warping). The additional correc- 
tional terms arise from the spatial dependence of u and A. 
Their derivation and the following analysis represent the 
subject and result of this paper. 

The effective Hamiltonian in Eq. Q can be simplified 
when Xp ^> I. Terms with \p\ 2 ~ kf can be dropped, 
given the length scales in this regime and the de Broglie 
relation. Now we wish to compare terms containing the 
potential u and gap A. To do this, we follow a simplified 
scheme to that defined in Ref. [17 , modelling the bilayer 
on a substrate as a parallel plate capacitor. 

Each layer of graphene has surface area A, and we take 
the dielectric constants of the material between the back 
gate and layer 1, and the bilayer, to be unity. Layer 1 
has charge Q = —riieA, while layer 2 has charge Q' = 
—ri2eA, where ri\ (712) is the density on layer 1 (2) (and 
n = ni + 712). The back gate and layer 1 are separated 
by a distance Lf,, while the two layers are separated by a 
distance cq. Applying a Gaussian surface around layer 1, 
the magnitude of the electric field is E = Q/e^A, where 
£o is the permittivity of free space. The voltage due to 
this electric field is thus QL^/eqA. The electric poten- 
tial energy due to the back gate (thus, the potential u) 
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Matching the wavef unctions at either side of the junction 
(■01 (0) = ^2 (0)), we obtain the boundary condition 
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Figure 3: Angular dependence of quasiparticle transmission through 
an n-p junction. 



is u = eQLb/eoA — n 1 e 2 Lb/eo. We assume that the elec- 
tric field from the back gate is screened poorly by layer 
1, so applying the same analysis to layer 2, we find that 
the magnitude of the electric field is E' = Q' /e A. The 
voltage produced by that electric field is V = E'cq, so 
the electric potential energy between the graphene layers 
(i.e., the gap) is A = n 2 e 2 c /eo. If we assume that the 
charge density is evenly distributed between the layers, 
n\ — n 2 = n/2, then u/A = L^/cq. With c$ ~ 0.3nm 
and Lb ~ 300nm, we find that u ^> A. By writing H c s 
in the form H c k = I A + a x B + a y C + cr z D, we compare 
each term and keep only the largest one in each group 
A, B, C, D. This produces an approximate Hamiltonian, 
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where rj £ {0, 1} and highlights the correctional terms. 

An n-p junction can be formed with two back gates, 
schematically shown in Fig. [3] Each gate can indepen- 
dently create an electrostatic potential over that region 
of bilayer graphene. Given our chosen length scales in Eq. 
([2]), we model the n-p junction as a Heaviside step function 
Q(x) — (1/2), with its derivative the Dirac delta function. 
Thus, u w (k 2 F /2m)[Q(x) — (1/2)], which also determines 
all additional terms in Eq. 

We define the problem in terms of plane waves on the 
left-hand and right-hand sides of the junction, ipi and ip2 
respectively: 



Ak x x 



(9) 



The Hamiltonian in Eq. ^ has plane and evanescent 
wave solutions, and the quasiparticles are chiral, such that 
when they pass from the conductance band at the left of 
the interface to the valence band at the right, k x changes 
sign [3] (see Fig. [2]). 

With the step defined to be at x — 0, we integrate Eq. 
m\ across it, Jq^{£ — H app ) dx and take the limit 6 — > 0. 



where the Fermi momentum kp = J k 2 + k 2 and 
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Using these equations, where k x = A— k 2 + 2m[(u/2) — e], 



= J-k 2 y + 2m[(u/2) 



e], K = Jk 2 + 2m[(u/2)-e] 



and k' = -Jkl + 2m[(u/2) + s], we calculate the trans- 
mission probability for a symmetric junction T(k y ) = |<i| 2 . 
We assume a wide strip, such that k y is invariant. We 
also set e — in the middle of the barrier for simplicity. 
Using k y = kp sin(0) (see Fig. [3]) we first calculate the 
transmission with only the leading-order terms in Eq. ( 10 1 



(by setting 77 = 0), finding agreement with Ref. [3] in that 
T(8) = sin 2 (20). Including the correctional terms from Eq. 
|l0| ) by setting 77 = 1, we obtain a correction to the inci- 
dent angle at which perfect transmission is seen (see Fig. 

Taylor expanding the full analytical result for T(9) 
around 77, we find that only the first-order term is impor- 
tant and obtain a potential-dependent result (providing a 
good fit up to u s» 50meV): 

2 v 

T(9) sin 2 (26») sin(46>) sin(0). (12) 

37i 

Assuming a wide graphene sheet (that is, a width w 
much greater than the length) and coherent quasiparti- 
cles, one can calculate the conductance from the trans- 
mission probability using the Landauer-Biittiker approach 
[T9] (taking into account two valleys and two spins), 
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With k y = 21X71 1 w (where 71 is an integer), we can write 
this as an integral and calculate it using the full numerical 
transmission probability, 





Figure 4: The dashed line shows the transmission probability without 
correctional terms applied. The solid line includes correctional terms. 
The dashed line shows perfect transmission at an angle of 45° to the 
interface, (a) Transmission at u = 40meV. (b) Transmission at u = 
80meV. Plotted with e = 0, V « c/300, 71 = 0.4eV, m = 0.035m e , 
A = 0. 
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for u — 40meV. This is a slight reduction from 2.12c 2 'wkp /nh 
for the case rj = 0. One can also calculate the Fano fac- 
tor [201 EI] (the ratio of shot noise to Poisson noise; for a 
review see Ref. jUj) numerically: 



jZj% d0cos(fl)r(fl)(i-r(0)) 

i:i 2 /2 decos(9)T(9) 



F 



0.241 



(15) 



for u = 40meV, showing a small increase compared to 
F = 0.238 when rj = (see Fig. [5}. 

In conclusion, we have extended the earlier derived 
low-energy effective Hamiltonian for bilayer graphene to 
incorporate a spatially dependent electrostatic potential 
consistently. We calculate the angle-dependent transmis- 
sion through an n-p junction and find T(6) = sin 2 (2f?) — 
(2u/37i) sin(40) sin(6*). Perfect transmission is still seen, 
but at a slightly increased angle. The conductance is slightly 
reduced to G = 2.1e 2 wkp/nh, whereas the Fano factor is 
slightly increased to F = 0.241 (both for u = 40meV). 
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